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Abstract. RenormaUzed versions of cosmological perturbation theory have been very 
successful in recent years in describing the evolution of structure formation in the 
weakly non-linear regime. The concept of multi-point propagators has been introduced 
as a tool to quantify the relation between the initial matter distribution and the final 
one and to push the validity of the approaches to smaller scales. We generalize the 
n-point propagators that have been considered until now to include a new class of 
multi-point propagators that are relevant in the framework of the renormalization 
group formalism. The large-A: results obtained for this general class of multi-point 
propagators match the results obtained earlier both in the case of Gaussian and non- 
Gaussian initial conditions. We discuss how the large-A: results can be used to improve 
on the accuracy of the calculations of the power spectrum and bispectrum in the 
presence of initial non-Gaussianities. 
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1. Introduction 

The understanding of tlie evolution of structures in tlie universe has improved 
dramatically with the emergence of cosmological perturbation theory in the beginning 
of the 1980 's [H El [3l H] and the possibility of performing high resolution N-body 
simulations [5l|6l[71|8]. The two approaches are complementary with perturbation theory 
being best suited for large scale calculations and N-body simulations working well on 
smaller scales. On intermediate scales the two should be compared as a consistency 
check on the approximations used in both approaches. The ultimate goal is to create 
computational approaches that can be used for sampling the cosmological parameter 
space and comparing with data from future galaxy surveys |9j and weak-lensing surveys 
This requires efficient numerical solutions and the computational cost of running 
N-body simulations makes it necessary to investigate semi-analytical approaches that 
can push the validity of perturbation theory further into the non-linear regime at small 
scales. 

The most relevant scales for parameter estimation are those where the baryon 
acoustic oscillations (BAO) are most prominent at wavenumbers of ~ 0.05 — 
0.25 h Mpc~^. At high redshifts standard perturbation theory describes these scales 
very well, but at the present time the one-loop results fail at A; ~ 0.1 h Mpc~^ and 
computation of the higher order corrections quickly becomes too time consuming (see 
[12] for two-loop results and jT3] for a thorough review of cosmological perturbation 
theory). New approaches are needed that either improve directly on the existing 
perturbation theory or formulate alternative solutions. 

The similarities between cosmological perturbation theory and quantum field theory 
have in recent years been exploited to create a diagrammatical interpretation of the 
equations analogous to the Feynman diagrams of particle physics. The formalism was 
introduced in [T3] and has been developed in different directions in [T5 | [T6 | ITT t ITS] . Other 
approaches that do not rely as heavily on understanding the field theoretical background 
that inspired them can be found in [T9l[20| [2T]. The general idea of these developments in 
the field is to rearrange the perturbative series and do a partial resummation of a certain 
class of diagrams to all orders. The validity of some of the approaches referenced here 
has been tested and compared against each other in |22] and against N-body simulations 
in [12]. While most approaches perform better than the standard one-loop result we 
are still interested in pushing the validity further into the non-linear regime. 

A key object in the diagrammatical formalism is the two-point propagator that 
describe the evolution of a Fourier mode over time. At large A;-values, i.e. small scales 
where standard perturbation theory breaks down, the perturbative corrections to the 
two-point propagator can be resummed to gain an analytical result that is well behaved 
on all scales. This was first achieved in |23] where the result was used to define a 
renormalized perturbation theory and later reproduced using the renormalization group 
(RG) formalism in pjj. In |2l] the RG approach was generalized to deal with non- 
Gaussian initial conditions and a large-fc result was found that included a non-zero 
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initial bispectrum. 

Multi-point propagators describe a combination of time evolution and couplings 
between Fourier modes. The n-point propagators were introduced in j25] in the context 
of renormalized perturbation theory and it was shown that they have a similar large-/c 
behaviour as the two-point propagator. In |26] the large-/c behaviour was obtained also 
in the presence of initial non-Gaussianities. In this paper we follow the RG approach 
of [T5| 121] to reproduce the large-fc results of [2^ with both Gaussian and non- 
Gaussian initial conditions. We also generalize the concept of n-point propagators to 
include a new class of multi-point propagators that is relevant in the RG formalism and 
obtain large-fc results for these. 

In Section [2] we present the equations that govern the evolution of the cosmic fields 
and define the multi-point propagators in the general form. Section [3] contains a recap of 
the RG formalism with non-Gaussian initial conditions and the recipe for constructing 
RG equations. In Section H] we compute the large-/c limit of the generalized multi-point 
propagators and in Section |5] we conclude and comment on our results. 

2. Evolution of Cosmic Fields 

We wish to follow the evolution of matter perturbations from some initial distribution 
in the early universe. We use the density perturbations 5 (x, r) defined through 

p(x,r) = p(r)(l + 5(x,r)) (1) 

and the peculiar velocity field v (x, r) of the fluid. These fields satisfy the continuity 
and Euler equations 

1^ + V ■ 1(1 + ^)v] = 
— + "Hv + (v ■ V)v = -V^ 

OT 

while the gravitational potential is sourced through the Poisson equation on subhorizon 
scales 

V^m = hl,,,{r)n\T)5{^,r) (3) 

Here the conformal Hubble parameter is defined as H = ^^jp^ = aH . We will assume 
that the background cosmology is Einstein-de Sitter with VL^ = 1 and Q\ = 0. 

In the regime where our approach will be valid, i.e., large scales, we can assume that 
the velocity field is irrotational so that the velocity divergence, ^ (x, r) = V ■ v(x, r), 
completely determines the velocity field. Taking the divergence of the Euler equation 
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and going to Fourier space we obtain 
9(5(k,r 



dr 



+ ^(k, t) + j d^^qd^p 5^(k - q - p)a(q, p)e(q, r)5(p, r) = 



^^ + 7/0(k,r) + ^7^25(k,r) (4) 
+ j d\d^v 5D(k - q - p)/3(q, p)^(q, r)0(p, r) = 
with the two momentum factors given by 

( . (q + p)-q ^ n, N (q + p)^q-p 

a q,P = 2 and /3 q, p = —- (5) 

These factors carry the non-hnearity of the equations by couphng different momentum 
modes of the density and velocity perturbations. 

We proceed following the approach of [18] and define a two-component field 
including the density perturbations and the velocity divergence 

^"<''''')^(::i^;:i)--'(x;;/«) 

where the new time parameter is the logarithm of the scale factor rj = ln(a/aj) and Oj 
corresponds to some initial time when the evolution was still in the linear regime on all 
relevant scales. In an Einstein-de Sitter cosmology the new time variable is equal to 
the logarithm of the linear growth factor, r] = ln{D/Di). It has been argued in |23] that 
this redefinition of rj can also be used in more general cosmologies such as ACDM even 
though it is not exact. For a scale independent linear growth factor the discrepancies 
affect only the decaying mode producing very small corrections to the end results at the 
relevant scales. 

With the field (pa the equations (jlj) can be combined to yield a more compact 
equation of motion for the fluid where repeated indices (running from 1 to 2) and 
momenta are being summed and integrated over respectively 

(Sabdri + ^ab) 0fe(k, T]) = e'^-fabci'k, -p, -q)06(p, v)^^^ v) (7) 

with 

I. ) (8) 



and the non-zero entries in the vertex factor are given by 

7i2i(k, p, q) = 7ii2(k, q, p) = |(5D(k + p + q)a(p, q) 
7222 (k, p, q) = (5D(k + p + q)/3(p, q) 



(9) 



The vertex describes how two Fourier modes will interact and create structure on a 
different scale in the process. 
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In the linear regime the right hand side of equation ([7j) is put to zero as there is no 
interaction between the different modes and we can solve for the linear propagator of 
the system gab-, that is the operator that evolves the linear solution 0° forward in time 
or the Green's function of the linearized system 



^a(k, r]a) = 9ab{Va, Vb)(i)biK Vb) for ria > r]b 



(10) 



By plugging this expression into equation ([7j) with jabc = we see that the linear 
propagator satisfies the equation 



{Sabd^, + Qab) gbciVa, Vb) = ^ac^oiVa " Vb) 

which can be solved to give 



gabiVa,Vb) 



3 2 
3 2 



+ 



2 -2 
-3 3 



Vb) 



(12) 



The last matrix term with the decaying exponential factor controls the linearly decaying 
mode while the first term controls the linearly growing mode. By selecting appropriate 
initial conditions we can follow the evolution of each mode choosing (pa proportional to 
Ua = (i) for the growing mode and Va = (-3/2) foi' the decaying mode. Due to the 
factor of e~'' in the definition of the two-component field in equation the first term 
in the linear propagator reduces to a time independent constant matrix in contrast to 
the growing mode of the linear propagator of [23] that grows as e^. 

2.1. Path Integral Formulation 

In the path integral formulation the key object is the generating functional Z which can 
be used to generate correlation functions between the fields. It is necessary to introduce 
an auxiliary field Xa{^,v) that will couple to the initial conditions (see [18] and |24] 
for a thorough introduction). In the case of Gaussian initial conditions the generating 
functional is 

Z[Ja,Kb,P^] = j VcpaVXb exp 1^ - i j dr]adr]b XaPab^{Va)S{Vb)Xb 
+ i J dr] [(paQabXb - e'^labcXaMc + Ja(t)a + K^Xb] 

where P^^ is the initial power spectrum that completely specifies the initial Gaussian 
matter distribution and Ja and Kb are sources for the two fields 0^ and Xb respectively. 

If we go back to the linear limit (setting e'^'jabc = 0) the integration over the two 
fields can be performed explicitly yielding 



(13) 



Zo[Ja,Kb;P^] =exp ( - / drjadrjb 

+ iJa(k, ria)gab{'na, 'nb)Kbi-k, r]h 



-JaiK Va)Pabik; Va, Vb)Jb{-K Vb) 



(14) 
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where we see that differentiating twice with respect to J gives the hnearly evolved 
power spectrum, P^i,{k; r]a, r]b) = gac{r]a,0)gbd{Vb,0)PcdW, and differentiating once with 
respect to J and once with respect to K gives the hnear propagator gab- These two 
objects along with the vertex 'jabc constitute the fundamental parts of the Feynman 
diagrams we will consider. They will be represented by the diagrams in figure [H where 
solid lines represent and dashed lines represent x- 

In Pl] it has been shown that in the presence of initial non-Gaussianities equations 
([Hj) and ( 1T3|) take the same form with equation (IT^ replaced by 



where 5^^^ is the linearly evolved bispectrum, [T^] represents a similar term with the 
linearly evolved trispectrum and the dots represent terms from higher order statistics. 
Equation (fT5|) shows that differentiating Zq n times with respect to the source will 
give 



Zo[Ja,Kb]P^,B^,- ■ ■] =exp|^-^ j dr]adr]b JaiK'na)Pabik;'na,'nb)M-^^'nb) 

~\ j (^Vadr]bdr]c 5^fec(ki, ^2, ks; r]a, rjb, r]c)Ja{'ki, r]a)Jb{'k2, r^^) Jc(k3, r?c) (15) 
+ [T^] H i J dr]adr]b Ja{k,r]a)gab{Va,Vb)Kb{-'k,r]b) 







Zo 6Ji{ki,r]i) ■ ■ ■6Jn(kn,Vn) 



= Sd Sa.-aA'^ir ■ ■ ,K;Vl,--- ,Vn) 

(16) 



where S^^...^^ represents the ra'th order statistic linearly evolved in time. The linearly 
evolved statistics will constitute additional fundamental building blocks in the Feynman 
diagrams and will be represented by a box with the appropriate number of legs as shown 
in figure [1] 

In equation fll3p we can replace the first term in the exponential function by the 
infinite series 

—XaiK 0)P0 (^)X6(k, 0) - ^5l(ki, k2, k3)Xa(kl, 0)Xfe(k2, 0)Xc(k3, 0) 

+— Tl,(ki, k2, ka, k4)xa(ki, 0)x6(k2, 0)xc(k3, 0)xd(k4, 0) + ■ ■ ■ 

including the initial bispectrum B^f^^, trispectrum Tj^bcd and higher order statistics, all 
coupled to the x-field. The fact that the initial statistics are coupled directly to the 
X-field only tells us that the auxiliary field plays a crucial role in connecting the evolved 
(/)-fields, (pirj), to the initial conditions. 

Having defined the fundamental building blocks of the Feynman series we can in 
principle build all the non-linear results by taking loop corrections into account order 
by order. In this work however we will follow the renormalization group (RG) approach 
of [18j which sets up a differential equation for the full non-linear objects that include 
a resummation of a series of diagrams to all orders in perturbation theory. 
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-D- 



Vb 




power spectrum: P^f^{k]ria,rih) 



bispectrum: B^^^{ka, kf,, k^; rja, rjh, r/J 



Figure 1. Fundamental building blocks of Feynman diagrams. The dashed lines 
represent the x-field while the full lines represent the </)-field. Time flows right to left 
in the propagator and the dots represent the series of higher order statistics. 



2.2. Generating Functional 

We begin by defining the generating functional of connected diagrams W 

W[J,,K,] = -tlog Z[J,,K,] 



[18) 



and tfie generating functional of one-particle-irreducible (IPI) diagrams F, which is 
given by the Legendre transform of W 



r[0a, Xb] = W[Ja, K,] - j dr/d^k ( J„0„ + K,Xb) 



(19) 



In this equation cf) and x actually represent the expectation values of the fields, but to 
keep the notation simple we will suppress this detail. We only note that the physical 
situation with Ja = Kf, = corresponds to the expectation values (0a) = (xa) = 
because we are dealing with fluctuations away from the average of the density and 
velocity fields. In the presence of the sources the expectation values of the fields are 
given by 



6Jn 



Xb[Ja,Kb 



5Ku 



(20) 



The IPI diagrams correspond to full vertices with any number of legs coming in 
from the right and going out to the left corresponding to differentiations of F with 
respect to (p and x respectively||| The derivatives are always taken at 0^ = Xa = to 
recover the physical situation. As discussed in |T8j one can realize that derivatives of 

I Notice the difference in notation as compared to and [5^ where F represents the multi-point 
propagators. 
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r with respect to (pa only, i.e. diagrams where every external leg represents a (/)-field, 
will vanish due to a closed loop of propagators. The second derivatives of F can then 
be written as 

^b) = 9baiVa, Tib) - ^^aXbi^] Va, Vb) ^^1) 
^S-/-.^^' ^b) = QabiVa, Vb) - ^Xa't'tik] Va, Vb) 

^?lxb(^'^^'^^^b) = iPab{k)S{ria)6{rib) + i<l>ab{k;Va,Vb) 

where the linear part can be read off of equation (IT^ and the $ and S terms are the 
self-energies arising from non-linear interactions. The notation F^^^^ represents the 
derivatives taken at 0^ = x;, = and includes a delta function in the momenta 



(2) x_ 5 r[0„,x6] 



Sn{k + k')rl\{k;Va,Vb) 



6(j)a{k, 1]a)6(f)b{'k',1]b) 



(22) 

'l>a,Xb = 



Along the same lines as in equation (12T]) the full power spectrum and propagator 
can be defined in terms of the second derivatives of W with respect to the sources J 
and K 

W^'Xik, r/,; k', r]b) = ^6n{'k + k')Pab{k- r/,, r/,) 

^'^ ^b) = -Snik + k')Gab{k; Va, Vb) ^^3^ 
Va; k',rib) = -Snik + k')Gbaik; rja, Vb) 
Wl^l^^ik,r^a;k\vb) = 

(2) 

The notation Wj j represents the derivatives taken at = = 0. The component 

(2) 

^KaKt being zero is a consequence of the fact that the two sets of second derivatives in 
equations (12T]) and f l23|) constitute inverse matrices of each other. This can also be used 
to express the power spectrum as a sum of two contributions 

Pab = Plb + P^b (24) 

given by 

Plbik-, Va, Vb) = Gac{k; Va, 0)Gbd{k; Vb, 0)P^,ik) 

pb (25) 

P^i,{k;r]a,Vb) = dsi ds2Gac{k;r]a, Si)Gbd{k;r]b, S2)^cd{k; Si, S2) 
Jo Jo 

where P^ is the equivalent of the linearly evolved power spectrum using the full 
propagator, and P^^ represents a new contribution from the interactions between the 
Fourier modes. It is clear from equation (125|) that the full propagator and the self- 
energy $ are key objects in this approach. Diagrammatically we will represent the full 
propagator with thick lines and $ with a circle with two external x^legs. 
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The full propagator, power spectrum and IPI n-point functions appearing in 
equations <^T^ and ( 12^ represent a generalization of the corresponding linear quantities 
in the sense that they include all loop corrections to the given quantity. A Feynman 
diagram drawn with the full quantities will represent an infinite sum of diagrams in 
standard perturbation theory. The resummed quantities have been exploited in [15j and 
[23j to build a renormalized perturbation theory of cosmological structure formation. 

The full bispectrum, trispectrum and higher order statistics can be obtained by 
taking higher order derivatives of W with respect to J, differentiating thrice for the 
bispectrum, four times for the trispectrum and so on. 



2.3. Vertices 

Similar to equation (12T]) we can identify the full vertex with derivatives of F 

= - 2(5(si - s)5{s2 - s)5{s-i - s)e"7a6c(k, -q, -k + q) (26) 
+ loop corrections 

The factor of 2 comes from the cfP' in equation (IT^ and the delta functions account for 
the fact that the tree level vertex is instantaneous, i.e. it has no time dependence. The 
full vertex will be represented by a filled circle in the Feynman diagrams so that the 
diagrammatic version of equation f l26|) reads 



2 \^ ^ (' ^ perm.)) 

+ (two-loop diagrams) + • • • ^^^-^ 

where I have suppressed the delta functions in the first term. 

As mentioned in Section 12.21 this definition can be generalized to higher order 
vertices with any number of legs coming in from the right and going out to the left 




(28) 

with the constraint that the sum over momenta on each side must be equal and opposite, 
k = — XliLi = 'YldLi Pi- "^^^ vertices defined in this way express how a number 
of Fourier modes can interact and combine to create power on different scales than 
before. This is entirely a non-linear effect as the linear evolution will preserve the 
initial statistics. 

Vertices with just one outgoing leg on the left side are the basis of the n-point 
propagators introduced by Their tree level expressions can be constructed using 
only the tree level vertex and the linear propagator, i.e. with no couplings to the initial 
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Figure 2. Decomposing a generalized vertex into n-point propagators. The red dotted 
lines indicate that the initial statistics work as couplings between individual vertices 
with just one outgoing x^leg as described in the text. 



conditions. On the other hand constructing diagrams with more than one outgoing leg 
requires the use of this kind of couphngs. By considering the possible diagrams one 
can realize that in fact they are constructed by gluing together a number of vertices 
with just one outgoing leg by coupling them to the initial statistics. See figure [2] for an 
example. 

The vertices constructed in this way are not fundamental objects in the 
renormalized perturbation theory of [22], but in the RG approach we will see that 
all the possible IPI vertices will in principle play a role, so it is relevant to gain some 
insight into the general behaviour of these objects as well. Before turning to the RG 
approach we will clarify the connection between the formalism and notation presented 
here and that of ESI and 1261. 



2.4- From Vertices to Multi-Point Propagators 

To get from the IPI diagrams in Section I2T^ to the multi-point propagators equivalent 
to those in |25] we need to attach propagators to all the legs of the vertices so that the 
resulting diagrams describe how the initial Fourier modes at a time r/b evolve and affect 
each combination of Fourier modes at a later time rja 

Va Vb 



(29) 

We will denote the multi-point propagators as with n being the number of 

incoming legs and m the number of outgoing legs as shown in the diagram above. As 
an 
be 



(2 1) 

an example the formal expression for V"^^^ in terms of full propagators and vertices will 



(ka,?7a; k5,kc,?7fe) = / Asx^S2^S^Gad{ka\Va,Sx 
J Vh 



% (30) 

'^xl^e<t>f(^<^^ ^b, S2\ kc, S^)Geb{h\ S2, Vb)Gfc{kc\ S3, Vb) 

with the implicit requirement that k^ = — k;, — kg. In this work we will distinguish 
between n-point propagators that have m = 1 and multi-point propagators that have 
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m > 1. This is in contrast to [25j and [26] where they only deal with the m = 1 case 
and the two terms are used interchangeably. 

(2 1) 

The n-point propagator V^^^'^ in this work is the equivalent of the three-point 
propagator F^j,^ of [25] . Apart from the differences that arise due to the extra factor of 
e~'^ in equation ([6]) compared to the two-component field of |SJ, the two propagators 
describe the same overall evolution. Naturally this equivalence holds between all the n- 
point propagators F*^""^^ of and the corresponding ^(""^'i) of this work. It was shown 
in [25] for Gaussian initial conditions and in [26] for general initial conditions that the 
n-point propagators are connected to the integration kernels of standard perturbation 
theory in such a way that the observable statistics at later times can be constructed 
directly from the n-point propagators. This is done in the same way as the vertices in 
Section [2l3] with more than one outgoing leg are constructed by gluing together n-point 
propagators via couplings to the initial statistics. A few of the contributions to the 
power spectrum will be 

(31) 

where the red dotted lines again indicate where the diagrams can be split into separate 
n-point propagators. It is clear from this direct approach where all the observables 
can be constructed from the ra-point propagators, that there is no need to consider the 
multi-point propagators with m > 1 in their framework. As described in Section 12.31 
this is not the case in our resummation scheme, so we will also consider the new class 
of multi-point propag ators with m > 1. 

We will conclude this section by studying the tree level expression for V^j^^' . The 
tree level version of equation fl30l) is obtained by replacing the full propagators with the 
linear ones and by use of the first term on the right hand side of equation (12^ for the 
vertex 

K?c,toee(k, Va] ki, ks, r];,) = -2 / ds e'gadiVa, s)7de/(k, ki, k2)^efe(s, r]b)9fcis, r/fe) (32) 

Vb 

The factor of 2 expresses the two possible ways of contracting the vertex with the two 
last propagators and the overall minus sign is due to the sign convention adopted in 
equation f|T3|) . 

To compare with [25] we will compute the component with a = 1 and growing mode 
initial conditions 

^lUtree(k, Va, ^1, ^^2, Vb)UbUc = "2 




/5 kix 2x^\ 
\7^2h^2k'i^~) 
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where x = (ki ■ h2)/{kik2). We see that setting r/b to zero and muhiplying by a factor 
e''" will reproduce the tree level expression of apart from the overall factor of —2 as 
expected. 



3. Renormalization Group Formalism 



The renormalization group formalism presented in [18] consists of constructing 
differential equations for the final statistics and full propagators introduced in Sections 
12.21 and 12.41 that gradually include more and more loop corrections from the non- 
linearities. The resulting differential equations can in some cases be solved analytically 
yielding a resummation of certain loop diagrams to all orders. This is particularly useful 
in situations where standard perturbation theory does not produce convergent results 
at any finite order because the non-linear couplings are large. As can be seen in the 
explicit expression of a and f3 in equation ([5]) the couplings in cosmological perturbation 
theory can be very large when the matter density is high. This happens at small scales 
(large momenta) when the initial distribution of matter has had time to collaps onto 
high density peaks. 



3.1. RG Equations 

To generate the renormalization group differential equations we need to introduce a 
cut-off that allows only the low-order loop corrections to contribute. This requirement 
translates into allowing low fc-values to contribute and can be realized by modifying the 
initial conditions with a step function as in [TB] 

P\k) ^ P°(A;) = P'ik)eH{\ - k) 
E°(ki, k2, ka) ^ 5°(ki, k2, ks) = 5°(ki, k2, k3)^^^,(A - {k^ + k2 + k^)) (34) 

The step function, Oh-, equals unity for k < X and zero for k > X. For the higher order 
statistics there are no unique way to introduce the cut-off. Here I have used the sum 
of the magnitude of the wavenumbers to define the cut-off scale. Another possibility 
could be to use maxjfcj}. The choice of cut-off scale will affect the behaviour at low 
A-values, but as we take the limit A — ?■ oo there should be no difference between the 
choices. 

Replacing the initial statistics with the A-dependent ones in the generating 
functionals and differentiating with respect to A will create the differential equations 
we are interested in. These can then be integrated from A = to A — )■ oo to regain 
the full dynamics encoded in the initial conditions. We will denote the A-dependent 
quantities by a lower case A as in equation flMj) . 

Because the cut-off appears only in connection with the initial statistics it is easy 
to calculate the A-derivative of Zx = Z[Ja, Kb, P^, B^, ■ ■ ■]. Starting from the equations 
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V(j)aVxb exp(- • ■ ) - - J dr]a,bd qS{r]a)5{7]b)xaXbS{X - q)Pab{(l) 

d?7aAcd^qi,2,3 ^{va)S{rib)6{r]c)xaXbXcS{x - E?i)^afec(qi'q2>q3) 



+ {r°} 



(35) 



1 



- j dVa,bd\SiX-q)P^b{q)S{Va)S{Vb)^^^^^^ 

+ {T°} + ■ ■ • 



where we have introduced the compact notation drja^b = (^Va^^Vb and d^q^ 23 = 
d^qid^q2d^q3 and the argument of the exponential function in the first line can be 
read off of equations (fT3|) and ( 1T7|) . 

The RG equation for W can be obtained by taking the A-derivative of equation 
(fT8|) . dxWx = —i-^d\Z\ and combining with equation fl35|) . The functional derivatives 
of Zx with respect to the source field K can be translated into functional derivatives of 
W\ by use of the identities 



5^Wx 
6KJKb 
6^Wx 
SKJKbSK, 



+ i 



1 6^Zx 
ZxSKJKb 
. 1 S^Zx 
'ZxSKJKbdK, 

A 



1 6Zx 1 SZ, 



A 



S^Wx 



Zx6KaZx6Kb 

+ XaXbXc 
A 



. 1 5'Zx 
'Zx6Ka6Kb 



IXaXb 



+ Xb 



5^Wx 



+ Xc 



SKaSKb 



(36) 



where we have used equation (120]) to rewrite single derivatives of W with respect to K 
in terms of the x^fisld. The terms that appear on the right hand side of the identities 
after the first term will be every combination of ^-fields and functional derivatives of 
Wx of order 2 or higher that corresponds to the number of differentiations on the left 
hand side. The coefficient of the first term will always be —i while the coefficients of 
the following terms are determined by an overall factor of i multiplied by an i for each 
X field and each Wx appearing in the term. The end result for the RG equation for Wx 
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IS 

dxWx = -j dr,,,,d\ 5{\ - q)P'M)KVa)5{Vb) [^XaXb + 



S^Wx 5^Wx \ 5^Wx 



(37) 



A VV X 

X I - XaXbXc + i I XaTTTTTr + XbTTTTTT + XcTTTTTr I + 



+ {r°} + • • • 

The main points to be drawn from this expression is that exactly one initial spectrum 
appears in each term on the right hand side and that all the functional derivatives of Wx 
are with respect to K and not J. This last point ensures that the functional derivatives 
will not generate more couplings to the initial statistics than the ones already explicitly 
present in equation fl?r|) . 

From equation (fT9|) we see that dxTx = dxWx, so the RG equation for Fx will be 
similar to equation (137|) and we can use that equation to generate RG equations for any 
object of interest in the theory. We have not yet set the sources to zero to recover the 
physical situation so we keep all the terms in the parentheses on the right hand side 
even though they will vanish in the physical limit. 

3.2. RG Equations for Multi-Point Propagators, Gaussian Initial Gonditions 

We start by considering the case of Gaussian initial conditions where only the first line 
of equation fl37p contributes. The two-point propagator Gab has already been studied 
using the RG formalism in fT8\. We will recap the results obtained there and then make 
the generalization to multi-point propagators. 

The RG equation for the two-point propagator is obtained through the definition 
of Gab in equation 



Ja,Kt=0 (38) 



o\Gab,x{k,'na,'nb) =- 



6Ja(k,r]a)5Kb{-'k,rib) 



2 



dr^edr/rfd^q 5(A - q)P'^iq)SiVc)5ivd)W^'^ 



where the sources have now been set to zero. Performing the differentiations on W the 
result can be written in terms of full vertices and propagators 

^j%,K,K„xi^^ Va] - k, r]b; q, Vc] -q, Vd) 

= j dsi ■ ■ ■ ds4 Gae,\{k] r]a, Si) Gfb,x{k; S2, r]b)Ggc,\{q] S3, r]c)Ghd,\{q; ^4, Vd) 

u) f (39) 

^xUf't's'Ph,)^^'^' -^1' -^2; q> -53; -q, ^4) - 2 / dsgdsg Ga^xik - q; s^, Sq) 

■ r£^<^,<^„A(k,Si; -q, S4; -k + q, S5)T^^l^^^.^^{k - q, Sg; q, S3; -k, S2) 
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The combination Ggc,x{q] s^,Q)Ghd,\{,(l] S4^,Gi)P^{q)5{\ — q) will arise naturally as 
the integration kernel in all RG equations, and diagrammatically it will be represented 
by a crossed square with two legs so that equation ( 138|) has the diagrammatical form 

A = 1 — h - — ^^^^ — 

dA 2 2 ^4Q) 

The above equation represents a general rule for constructing RG equations in the case 
of Gaussian initial conditions, where only one-loop diagrams are present on the right 
hand side. This simple recipe arises from the fact mentioned in Section 13.11 that only 
one initial power spectrum appears on the right hand side of equation (137|) and loops 
can only be constructed by use of couplings to the initial statistics. 

The general rules as set up by [18j for generating the right hand side of the RG 
equation for a given quantity in the Gaussian case are 

• Draw all one-loop corrections to the quantity using the full A-dependent 
propagators, power spectra and vertices. 

• Perform the A-differentiation of the full expressions by considering only the explicit 
A-dependence appearing in the theta function of P°. 

We can use these rules to generate the RG equation for the three-point propagator 
V^^j^ . By taking the A-derivative of equation (130|) we see the general structure of the 
RG equations we will encounter 

(2 1)/ / (3) (3) 

^>yabc,\ = J dSi,2,3 [dx iGad,\) ^ Xd'Pe'Pf ,\^^b,\G fc + Gad,X^ ^x^X {Geb,x) Gfc,X ^^^^ 

+ Gad,X^'^l^^^^^x'^eb,xdx {Gfc,x) + Gad,xdx {^xl'f>e<i^f,>) Geb,xGfc,x}j 

The first three terms on the right hand side represent the RG corrections to each leg 
as already expressed in the RG equation for the propagator in equation (HOjl . The last 
term represents the corrections that affect the vertex directly and using the above rules 
we can express the vertex RG equation diagrammatically as seen in figure [31 

This approach generalizes to higher order multi-point propagators in the sense 
that the RG equation can always be separated into one-loop diagrams that represent 
corrections to the two-point propagators in the legs and those that represent corrections 
to the vertex. 

The right hand side of the RG equations for vertices with more than one x^^g 
will in addition to the one-loop diagrams similar to those in figure [3] contain tree level 
diagrams of the form 

S*^o 

(42) 

Due to the delta functions 5{\ — q) from the integration kernel and 5(q + kj) from the 
propagator these diagrams will only contribute in a single point X = ki and will have no 
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d 




Figure 3. RG equation for the vertex F^" 



effect on tlie final solution of the differential equation. Thus we can safely ignore these 
deviations from the rules stated above. 

RG equations for the final power spectrum, bispectrum and higher order statistics 
can be generated in a similar way as in equation (130|) . For the power spectrum we see 
from equation (125|) that the only missing piece is the RG equation for which will 
include diagrams like 



where we see that vertices with more than one outgoing x^^g show up for the first time. 
To calculate the power spectrum we encounter vertices with up to two x-legs, for the 
bispectrum we would have vertices with up to three x^l^gs and so on for higher order 
statistics. This is why it is relevant to generalize the n-point propagators of [25j to 
multi-point propagators with more than one outgoing leg in the RG formalism. 

As a final note we see that we can obtain the one-loop result from ordinary 
perturbation theory by keeping all quantities in the diagrams at tree level and integrate 
the resulting expression over A. In this case only the first three diagrams on the right 
hand side of figure E] will contribute due to there being higher order vertices in the last 
four. 

3.3. RG Equations for Multi-Point Propagators, Non-Gaussian Initial Conditions 

With the presence of initial non-Gaussianities the RG equations become more 
complicated. The first line in equation (137|) will generate the same diagrams as presented 
in Section [3l2| but in addition to these the remaining terms will generate an infinite series 
of diagrams with up to two loops for the bispectrum, three for the trispectrum and so 
on. The integrations over will give rise to higher order integration kernels of the form 



G'.rf,A(gi;si,0)G,e,A(g2;s2,0)G,/,A(g3;s3,0)i?,%(qi,q2,q3)(5(A-Egi) ,[T1, ■■■ (44) 




and 




(43) 



We will represent these diagrammatically as a crossed box with the appropriate number 
of legs for each order of the initial statistics. In accordance with the comments after 
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-..c^^::?:^^.. _..<^.v.. 



Figure 4. All loop diagrams derived from the second term of equation (j37p for the 
two-point propagator. 



equation (137|) the general rule will be that only one coupling to the initial conditions in 
the form of the integration kernels will be present in each diagram. Thus the general 
rules from Section 13.21 for constructing the right hand side of the RG equation for a 
quantity can be reformulated in the non-Gaussian C3jSG clS 

• Draw all loop corrections to the quantity with exactly one crossed box using the 
full A-dependent propagators, statistics and vertices. 

The differentiation with respect to A is already taken care of by including the crossed 
box. Figure H] shows all the diagrams arising from the bispectrum integration kernel for 
the two-point propagator. As can be seen all diagrams are two-loop corrections and 
similarly the trispectrum will only generate three-loop diagrams and so on. This will 
also be the case for the RG equations for vertices with just one x^^g. 

For vertices with more than one x^l^g there will be additional diagrams at lower 
order arising from each integration kernel, but the delta functions in the wavenumbers 
will limit the range at which they contribute. For the vertex "^-^^^^^^ the bispectrum will 
generate six additional one-loop diagrams 



When performing two of the momentum integrations in equation (137|) we get q2 = — ki 
and qg = ki — q^^, where ki is the momentum associated with the lower leg to the left 
in the diagrams shown in equation (H5|) . The delta function in equation (jH]) will then 
reduce to 5{\ — gi — fci — |ki — q^|), so if ki is much larger than the scale at which 
the bispectrum has a significant amplitude, that is ki ^ gi, the diagrams will only 
contribute in a very short range 2ki < A < 2ki + 2gi. The same arguments apply to 
loop corrections from the higher order statistics, thus in the large-A; limit we can neglect 
the contribution from the diagrams at lower order in the perturbation series arising from 
each initial spectrum. 

It should be noted that the above simplification of the right hand side of the 
RG equations applies only in the large-fc limit. Ultimately we will be interested in 
calculating the final statistics at low fc-values where the lower order diagrams cannot 
be neglected. In particular the power spectrum will get a new one-loop contribution 




+ (1 perm.) + 




' » ^ 

+ (Iperm.) + ^^f^ 1" (Iperm.) 




(45) 
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from the initial bispectrum given by the last diagram in equation fl3T]) that should be 
included in the calculations. We will now turn to the large-fc limit of the multi-point 
propagators and neglect these lower order diagrams. 



4. Large— A; Limit of Multi— Point Propagators 

A first approximation for solving the RG equations would be to keep everything on the 
right hand side at its linear level. As mentioned in Section 13.21 this will reproduce the 
one-loop results from standard perturbation theory, so we will be interested in going 
beyond this level of approximations. 

The large-fc limit of the n-point propagators has been obtained both in 

the Gaussian |25] and the non-Gaussian case [26j . Here we reproduce their results using 
the RG formalism and show that a similar large-fc behavior appears for the new class 
of multi-point propagators V^("''") , m > 1 introduced in Section 12. 4[ 



4.1. Large-k Limit with Gaussian Initial Gonditions 

In the Gaussian case the large-A; regime is defined by kiCj^ ^ 1 where al is the initial 
velocity dispersion 



<yl = -J d^q (46) 



and ki are the momenta in each leg of the multi-point propagator. In this regime the 
integration kernel can be kept at linear order because the dominant g-values are low 
and taking the initial fields to be in the growing mode it reduces to 

El ~ gac{si, 0)gbd{s2, 0)UcUdP°{q)6{\ - q) = UaUb9n{si)9H{s2)P^{q)6{\ - q) 

(47) 

where the step functions ensure that Si,S2 > 0. We follow the approach of ^8j, where 
the large-A; limit of the two-point propagator has been obtained, and take the ordinary 
one-loop expressions for the right hand side of the RG equations as our starting point. 
Here we will go through a detailed calculation for the three-point propagator V'^'^'^\ 

4.I.I. Large-k limit of V^'^'^\ the Gaussian case Keeping the vertices at tree-level 
means that only the last diagram in equation (HOj) and the first three diagrams in figure 
[3] will contribute. This gives a total of 6 diagrams to consider, one for each propagator 
and three for the vertex itself. We will go through the calculations for the first vertex 
correction as an example. The complete diagram is 

ki , r]b 




k 



2 



(48) 
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Keeping everything at tree level and using equation fHTl) for the kernel we get 

-8 j dVsa.M e'"+"''+^"P°(g)5(A-g)e(s„)^(sfe)^ad(r?a,Sa)7de/(k,-q, -k + q)Me 

X gfg{sa, Sj,)7g/ii(k - q, ki + q, \^2)gic{s^, r]b) ^^^^ 

X ghj{sv, Sb)ijki{-'i<^i - q, q, ^i)ukgib{sb, vt) 

where the intermediate times are restricted by r]a > Sa > > Sb > rji, via the theta 
functions in the linear propagators. 

In the ki ^ q limit two of the vertices can be reduced further using the 
approximation 

k ' CI 

7de/(k, -q, -k + q)Me ^ ~^^df (50) 

and the linear propagators can be contracted with these delta functions so that the 
dependency on Sa and appears only in the exponential function in the first line of 
equation (H^ . Performing the two integrals over Sa and s^ then yields the large-/c result 

-2 j As^ e'^gag{ria,s^)-igM{\^MM)9ic{sv,rib)ghb{sv,rib){e'^^ - e"'')(e'" - e^") 

X U\^^P\q)5{X-q) ^^^^ 

J q q 

and the momentum integration in the last line can be rewritten as 

(k.k,)^Jd\^Si\-q) (52) 

Going through a similar calculation for the other five diagrams will change only the 
momentum factor in front of the integral in equation (l52l) and the limits of integration 
over Sa and Si,. Adding all six contributions gives 

^aV;£ A = j ds„ e^'-gagiva, s^)'yghiik, ki, k2)gicisy, r]b)ghb{sv, Vb) 

X (^e (e^^ - e^")' + {kl + kl + 2ki ■ ka) (e^" - e^'')' (53) 
- 2k ■ (kl + ka) (e^-^ - e^") (e^" - e"") )\ j ^\ ^-^^(^ ' l) 

Using the fact that k = — ki — k.2 we can reduce the sum in the parantheses to 
k"^ (e''" — e'"')^, and recognizing the tree level expression from equation (I32p the end 
result for the large-A; limit of the one-loop calculation of the RG equation for V^"^'^^ is 

9.V^S ^ -VlSL^^^^l j d'o ^*(A - (54, 

We can now go beyond the usual one-loop result by promoting the tree level expression 
for the multi-point propagator on the right hand side to the full A-dependent expression. 
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(3 1") fl 2) 

Figure 5. Tree level expression for V^^^"^^ and V"^^^ . 

The integration over A is then straightforward and in the hmit we are interested in, 
A — i- oo, the result is 

V^^^ = K.?i'Lexp (-^^ (e- - e")^) (55) 

where we have used the initial condition V^^l^^ = V^^^/tree '^v has been defined 
in equation (H6|) . This is the result obtained in [25] for the three-point propagator by 
resumming an infinite series of higher order loop diagrams in standard perturbation 
theory. The big advantage in the RG calculation is that we need only consider the 
explicit expressions of one-loop diagrams. 

The result shows that in the large-/c regime where standard perturbation theory 
is divergent at each order the full resummation of the dominant diagrams still gives a 
finite result with an exponential damping for large /c-values. 

4.1.2. Large-k limit of general multi-point propagators, the Gaussian case For the 
general multi-point propagators there are two additional considerations that should be 
noted. First the tree level expressions for the diagrams will be a sum over a few different 
realizations that should be studied separately. Secondly the propagators with more than 
one outgoing 0-leg will have couplings to the initial power spectrum as part of their tree 
level expression. See figure [S] for two examples. We will now show that a calculation 
similar to that just presented can be performed for each of these tree level diagrams. 

The lesson learned from the above calculation of one-loop diagrams is that in the 
large-/c regime each loop diagram can be reduced to the tree level expression multiplied 
by a momentum factor and a time factor as in equation (15T]) . The momentum factor 
can be rewritten as in equation fl52p where the dot product in front of the integral is 
determined by the two legs in the tree level diagram that the loop is attached to. The 
time factor arises from the integration over e*'""'"*'' where Sa and Si, are the two times of 
interaction where the loop hits the tree level diagram. The limits of the integrations 
can be read off from diagrams like those in figure |5] depending on which section of the 
tree level diagram the legs of the loop touch. 

Let us start by considering a segment of a tree level diagram with a power spectrum. 
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ki: 
ks: 

■h 

k4 



.-k5 



Figure 6. A tree level diagram for l/^^'^). 



There will be three types of loop diagrams with both legs attached to such a segment 



S2 



Si S2 







Sl S2 



Si 



(56) 



The momentum factor will be the same for all three diagrams while the time factors on 
the other hand are determined by integrations from the initial coupling to the power 
spectrum at time t = to either Si or S2 giving the sum 

^ (e- - 1)^ - (e- - 1) (e- - 1) + ^ (e- - if = ^ (e^ - e-)^ (57) 

Similarly any loop connected to the power spectrum segment with just one leg can be 
attached on either side of the box, but the sum of the two diagrams give a time factor of 
(e^2 — e^^) times the common time factor associated with the other leg and the common 
momentum factor. These considerations show that even though the segments with a 
power spectrum contribute with more loop diagrams than a propagator segment the 
extra contributions add up to the same result as a propagator segment will give. Thus 
we can treat any appearance of the power spectrum in the tree level expressions in the 
same way as a regular two-point propagator when we sum over all one-loop corrections 
on the right hand side of the RG equations. The two tree level diagrams for V^(^'^) in 
figure [5] can then be treated completely analogous to the V^"^'^^ calculation above. In 
general any tree level diagram for the multi-point propagator gous 
to a tree level diagram for so we will now restrict our attention to the n-point 

propagators. 

Each tree level diagram for the {n + l)-point propagator V^("'^) will contain n — 1 
mergings of legs on the right hand side of the diagram (see figure [H] for an example). 
The number of two-point propagators in such a diagram is 2n — 1 and the number of 
possible one-loop diagrams is n{2n — 1). We can treat all these one-loop corrections 
systematically by adding up diagrams with loops connected to a certain segment of the 
tree level diagram separately. 

If we consider a segment where two legs merge (as in the bottom part of the diagram 
in figure [6l where the legs with momenta k4 and ks merge between the times r]b and S4) 
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we can use the calculation for V^"^'^^ to add all the loop diagrams with both legs on the 
same segment. The resulting momentum and time factor will be 

|ki + kjf (e^' - e''")^ (58) 

where Si is the time at which the next merging takes place and kj and kj are the momenta 
before the merging. Similarly there are three loops with just one leg attached to the 
segment and the other leg on a propagator with momentum p going between times Sp^ 
and Sp2- These can be added to give the momentum and time factors 

2p • (ki (e^™ - e''") + kj- (e"'" - e"'') + (k^ + k^-) (e"' - e'^)) (e"^^ _ e'^^) 

= 2p • (k, + kj) {e" - e"") {e'''^ - e'"^ ) ^^^^ 

where we see that the time of merging, Sm, does not appear in the final result. This 
taken together with equation (158|) shows that the segment can be treated as just one 
two-point propagator with momentum kj + kj between the times and s; when we 
sum over one-loop diagrams. 

Using this approach iteratively for every merging of legs in the tree level expression 
under consideration we end up with just two legs merging at the final vertex where 
the result will depend only on the total momentum k = — ^ kj and the initial and 
final times rja and 77^. The intermediate times will only appear in the explicit tree level 
expression that will be common to all the one-loop diagrams. As argued above the 
approach works also for tree level diagrams with more than one outgoing leg on the 
left hand side where there will be one or more legs with a power spectrum. The only 
difference being a change of sign in the momenta considered so that k = ^ Pj = — ^ kj 
where Pj are the momenta on the left hand side and kj the momenta on the right hand 
side. 

When the sum over all one-loop corrections is taken the resulting momentum 
and time factors will be independent of which particular tree level realization we are 
considering so adding all possible one-loop diagrams for a given multi-point propagator 
will give 

^aK; '.aibi...6„,A = -k ^ 3 J d'q ^^^(A - 9) 2^ (tree level diagrams) 

(60) 

which is identical to equation (^^. Promoting the tree level expression on the right 
hand side to the full A-dependent one before integrating, we get the general large-/c 
result 

Tr{n,m) Tr(n,m) f ^ '^v ( rja 'Jii'l^^ ^P,^\ 

This result has been obtained for m = 1 and ?7{, = in ^25j. The fact that it generalizes 
to m > 1 is not obvious from the construction of this class of multi-point propagators. 
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As explained in Section 12.41 they can be constructed by combining two or more 
n-point propagators via couplings to the initial statistics, but the result is not just a 
product of the V^("~^'^) results. There are additional couplings between the individual 
n-point propagators that modify the result. 

Let us consider a specific diagram for to illustrate this point 



Va. Pi 



Vb 



Vb 



(62) 



Here we have a two-point propagator v}^'^^ 



ab 



Gab running from rjf, to rja with momentum 

(2 1) 

Pi, and a three-point propagator V^^j- running from rjt, to rja merging the two momenta 
Pi and — k into The large-fc results for these two have the damping factors 



exp 



and 



exp 



oVb} 



respectively while the combination in the diagram in equation 



is damped by 



exp 



(63) 



(64) 



where k = p^^ + pg. The additional couplings that arise in the combination of the two 
n-point propagators correspond to the term Ip^ ■ p^ in k'^ . 

This is analogous to bremsstrahlung processes in particle physics where couplings 
to an external potential leads to emissions or absorptions of particles. Here the initial 
conditions take the role of the external potential while each Fourier mode acts as an 
individual particle and the coupling to the initial conditions induces couplings between 
the two modes. 



J^..2. Large-k limit with Non-Gaussian Initial Gonditions 

The large-fc assumption of Section 14.11 is equivalent to assuming that the momenta k, 
running along the legs of the tree level diagrams are much greater than the momentum 
q in the power spectrum integration kernel defined in Section 13.21 We will also use this 
assumption in the case of initial non-Gaussianities so that the momenta kj are always 
much greater than the loop momenta q^- in the integration kernels of equation (l44l) . 

As seen in equation (l37j) the right hand side of the RG equations can be split into 
separate contributions from each order of the initial statistics. The contribution from 
the power spectrum will not change due to the presence of initial non-Gaussianities so 
we can reuse the results from Section 14.11 directly. If we do not rewrite the momentum 
factors as in equation (152|) the end result in equation (160|) can be expressed as 
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where we have promoted the tree level expression to the full A-dependent one. 

Next we turn to the contribution from the bispectrum. As we argued in Section 
13.31 we will only consider two-loop diagrams as the ones in figure HI Out of these only 
the three in the first line will contribute at tree level, and the first one will have three 
vertices coupling directly to the large momentum in the tree level diagram as opposed 
to just two for the last two diagrams. We can see from the vertex coefficients a and /3 
in equation ^ that these couplings introduce a hierarchy between the diagrams of the 
order 

k ■ qi k ■ k ■ q3 k ■ k ■ (q^ + qg) q2 ■ , . 

2 2 ~ ^ 2 \ \ r2 ~ (66) 

qi qi qi qi |q2 + q3r ^2 

so in the large-/c regime we need only consider the first kind of diagrams where all legs 
are connected directly to the tree level diagram. This argument applies equally well to 
the case of higher order statistics. 

Similar to equation (H7|) the bispectrum integration kernel reduces to 

^ Ma«bMc6'H(Sl)6'H(s2)6'H(s3)5°(qi, q2, f\^)5 {\ - Y. qi) , 

(67) 

when the initial fields are in the growing mode. This means that all three vertices that 
attach the bispectrum to the tree level diagram can be approximated using equation 
(150|) and the delta function in the indices ensures that we are again left with just the 
tree level expression times a momentum factor and a time factor as in equation fl^T]). 
The momentum factor will be 

/ d^q,d^q2d^q3 ^^^^5°(qi, q2, qs)^ (A - E^O (68) 

J q\ qi q-i 

where kj are the three momenta running along the segments of the tree level diagram 
where the loop is attached. 

The time factor consists of the integral J dsids2ds3 e^^"'"'^^^''^ where Si are the times 
of interaction with the loop momenta and the limits of integration depends on which 
segments the loops are attached to. 

These observations generalize readily to higher order statistics, so the task is now 
to consider how these momentum and time factors add when we sum over all diagrams 
on the right hand side of the RG equations. We start by considering a segment with the 
merging of just two legs in the tree level diagram which is equivalent to the three-point 
propagator V^'^'^^ . 

4-2.1. Large-k limit of V^^'^\ the non-Gaussian case For generality we will consider 
the nth order of the initial statistics. This means that there will be n loop legs that 
should be attached to the tree level diagram 



k 



Vb 



(69) 
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We will denote the momenta in the loop legs as ■ ■ ■ q„ and the times when they 
interact as Si ■ ■ ■ s„. It is not necessary to distinguish between the different loop legs 
since we integrate over both and Si and the initial statistics are symmetrized in their 
momentum dependencies. 

If there are m loop legs attached to a specific tree level leg in between times s^^ 
and the time factor from that leg is determined by the m integrals 



V2 r^m rS2 

dsm e^" / ds„_i e^™-^ ■ ■ ■ / dsi e'' (70) 



which by induction can be shown to yield 

dsm e^^T — ^—-T (e^'" - e^"i)^™-') = ^ (e'^^ - e^"i)™ (71) 
[m — 1)1 ml 

Thus if we have m loop legs attached to the left leg in equation f lB^ they will contribute 
with a time and momentum factor given by 



1 „ k ■ n 

m! ' -^-^ a 

1=1 ' 



(72) 



The remaining n — m legs should be attached to the two legs on the right hand side 
in equation (1^^ . The resulting n — m + 1 diagrams will have a common time factor 
of (e''" — e'"')""™', but the fraction connected to the time factor (as in equation (17T|) ) 
and the momentum factor will depend on exactly how many loop legs are connected to 
which tree level leg. With just one leg to attach we get two terms 

-ki ■ q, - ka ■ q, = k ■ q^ (73) 

where we have left out the factor of l/qf- With two legs one of the three terms can be 
rewritten using the invariance under q^ q^ so that 

^(ki ■ qi)(ki ■ q^) + (ki ■ qj(k2 ■ q^) + ^(ks • qi)(k2 ■ q^) 

= ^ ((ki ■ qi)(ki ■ q^) + (ki ■ qj(k2 ■ q^) + (ka ■ qi)(ki ■ q^) + (ka ■ qi)(k2 ■ q^)) 
= ^(k- qj(k- q^) 

(74) 

This result generalizes easily for n — m legs giving 



{n — m)\ ql 



1=171 



where we have reinstated in the factors of 1/qf and the common time factor. 

Equation (I75p taken together with equation 0721) gives the total momentum and 
time factor obtained by adding all diagrams with m loop legs attached to the left leg in 
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equation (169|) . Now adding all possible loop diagrams from the nth. order initial statistic 
results in a sum over these factors from m = to m = n 



m=0 ^ ^ j=l ^« /^gx 

1 " L- / \ 1 " L- 

= ^ n ^ E J (^''- - ^'■)"' - = f^"- - n ^ 

j=l ^« m=0 ^ ^ ■ j=l ^« 

where we see that just as in the Gaussian case we have eliminated the dependence on 
the intermediate time s„ so that it will only appear in the explicit expression for the 
tree level diagram in equation (169|) . 

The full result for the RG equation for V"*^^'^^ after promoting the tree level 
expression to the full non-linear one is then 



We can express the initial statistics in terms of correlation functions of the initial density 
perturbation 6o 

S{q, + q,)P'{q) = {5o{q,)So{q,)) 

^{^1 + q2 + q3)^°(qi> ^2, qs) = ('^o(qi)'^o(q2)<^o(q3)) (78) 



This enables us to rewrite the parentheses in equation fl77|) more compactly so that the 
end result after integrating and letting A — ?► oo is 

vS.» = KJiiL exp (t / d'q. ^*«(*))) (™) 

The change in sign for the power spectrum term compared to equation (177|) is due to 
the fact that the delta function S (q^ + q2) has not yet been enforced in this expression. 

4-2.2. Large-k limit of general multi-point propagators, the non-Gaussian case To 
generalize to all multi-point propagators we need to go through the same considerations 
as in Section UT] For the multi-point propagators with just one outgoing leg on the left 
hand side we can again use the calculation for V^"^'^^ above to show that every segment 
that merges two legs from the right can be treated as just one two-point propagator 
with the sum of the two initial momenta running along it. An iterative use of this gives 
the general result analogous to equation (179|) . 

On the other hand the case of multi-point propagators with more than one outgoing 
leg on the left hand side is more complicated in the presence of initial non-Gaussianities. 
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Figure 7. Example of a tree level diagram with a trispectrum merging 4 legs. The 
reasoning behind the last two diagrams is explained in the text. 



Apart from tree level segments that include an initial power spectrum there will also 
be segments that merge two or more legs by direct couplings to higher order initial 
statistics (see the first diagram in figure [7] for an example). We can split such diagrams 
into two parts by letting all the legs from the initial statistic go out to the time rjh before 
proceeding to the next interaction as seen in the second diagram of figure [71 Now any 
loop correction to such a tree level diagram that has a leg attached inside the center 
circle will not contribute to the RG equation when we sum over all diagrams, because 
the sum of the momenta coming from the initial statistic has to be zero. It is equivalent 
to the arguments leading up to equation (175]) except with k = 0. This means that we 
can treat the diagram as separate multi-point propagators all going from rjh to rja as 
shown in the last part of figure 13 

Each of these separate sections can again be reduced to a single two-point 
propagator carrying momentum and calculating loop diagrams along the same lines 
as was done for V^'^'^^ we arrive at the same final result as seen in equation (17^ 



V 



(n,m) 



ai---ambi---bn,\^oo 



(n,m) 



ai---ambi- 



,tree 



exp 



E 

,n=2 



, j=i 



(80) 



This is the result obtained in [26] for the n-point propagators V^'^~^'^\ that we have 
shown also holds for the generalization V"'-"'™^ with m > 1. 

The multi-point propagators can be interpreted as a measure of how well the initial 
conditions are preserved later in time. The linear two-point propagator for instance 
tells us that the linear evolution preserves the initial statistics perfectly because it is 
normalized to 1 at all times and scales. The full propagators on the other hand all 
show an exponential damping at large A;-values meaning that the memory of the initial 
statistics is erased at small scales where the non-linear collapse of the matter density 
peaks has proceeded far beyond the regime where standard perturbation theory is valid. 

The generalized results presented here show that even though the multi-point 
propagators m > 1 are in a sense more connected to the initial conditions 

through their tree level expressions, they lose memory of the statistics at the same rate 
as the n-point propagators. 
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It was noted in [26] that in an isotropic background universe the terms in the 
exponential of equation ( ISUj) cannot depend on the direction of k. We saw in the 
Gaussian case that the contribution from the power spectrum is proportional to k'^. 
Considering the term from the bispectrum the k-dependence is contained in the factor 

m = I d^qid^q,d^q3 ^^^^fi°(qi, q^, Qa) (81) 

We can perform an explicit integration over the orientation of k yielding 4:TTf{k) = 0, i.e. 
the bispectrum contribution vanishes. This will be the case for all the terms with uneven 
n in the sum. Only the even terms survive and the first correction to the Gaussian case 
comes from the trispectrum. 

5. Conclusions 

In this paper we have investigated the large-/c limit of multi-point propagators in the 
framework of the RG formalism presented in [18]. We have reproduced the results of 
|25] in the case of Gaussian initial conditions and [26] in the presence of initial non- 
Gaussianities for the special class of n-point propagators that was introduced in |25| . 
In addition to this we have obtained large-/c results for the generalized multi-point 
propag ators with m > 1 that play a role in the RG formalism along side the 

n-point propagators. 

All the multi-point propagators show the same exponential damping at large k- 
values that was first found for the two-point propagator in |23j for the Gaussian case. 
When initial non-Gaussianities are considered the first non-zero correction comes from 
the trispectrum which in most cases will only give minor deviations from the Gaussian 
result. In [26] a specific local model was considered with /nl = 10'^ and (^nl = 10^ that 
gave up to 2% weaker decay ai z = 0. With more realistic values for /nl and gf^L the 
deviations are much smaller. 

In [2Sj and [26] the large-fc limit of the n-point propagators is used directly to 
construct the power spectrum and bispectrum at late times and on small scales. The 
result is then combined with standard one-loop perturbation theory on large scales 
by interpolating between the two regimes. In the RG approach the large-fc limit of 
multi-point propagators could be used in a similar way as it has been for the two- 
point propagator in [18] where the time dependence from the large-A; result is assumed 
also to be valid for small fc-values. The large-/c results for the multi-point propagators 
represent a renormalization of the IPI vertices so an implementation of the RG equations 
that go beyond the tree level approximation for the vertices might be possible. 

Another possibility would be to calculate the power spectrum and bispectrum in 
the presence of initial non-Gaussianities. The framework to do this has been set up in 
Section HI but the RG equations become more complicated when we are also interested 
in the behaviour at small /c-values. In this regime the kinds of diagrams in equations 

and that could be neglected in Section H] play an important role. In particular 
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the power spectrum will get a new one-loop contribution from the initial bispectrum 
given by the diagram 

(82) 

It seems that including a diagram like this in the numerical solution of the RG equations 
should be possible. When it comes to the bispectrum the multi-point propagator 
provides us with the diagram 

':©-Bi— 

(83) 

that might work as a first approximation to the renormalized bispectrum. 

The approach followed in the calculations presented in this article corresponds to 
resumming only loop diagrams where all the loop legs are attached directly to the 
multi-point propagator in standard perturbation theory. In [27J it was shown that in 
the Gaussian case it is possible to resum another class of diagrams in the large-fc limit 
by essentially using a resummed power spectrum in the integration kernel of equation 
( H7|) . It would be interesting to see if their result can be reproduced in the case of 
multi-point propagators. 
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